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Based on the modified Landau-Lifshitz-Gilbert equation for an arbitrary Stoner particle under an 
external magnetic field and a spin-polarized electric current, differential equations for the optimal 
reversal trajectory, along which the magnetization reversal is the fastest one among all possible 
reversal routes, are obtained. We show that this is a Euler-Lagrange problem with constrains. 
The Euler equation of the optimal trajectory is useful in designing a magnetic field pulse and/or a 
polarized electric current pulse in magnetization reversal for two reasons. 1) It is straightforward to 
obtain the solution of the Euler equation, at least numerically, for a given magnetic nano-structure 
characterized by its magnetic anisotropy energy. 2) After obtaining the optimal reversal trajectory 
for a given magnetic nano-structure, finding a proper field/current pulse is an algebraic problem 
instead of the original nonlinear differential equation. 

PACS numbers: 75.60.Jk, 75.75.+a, 85.70.Ay 



The advent of miniaturization and fabrication of mag- 
netic particles with single magnetic domains [H, called 
Stoner particles, makes the Stoner- Wohlfarth (SW) 
problem [2J very relevant to nano-technologies and nano- 
sciences. One current topic in nanomagnetism is the 
controlled manipulation of magnetization dynamics of 
Stoner particles. Magnetization state can be manip- 
ulated by a mag netic fieldH i, H through the usual 
magnetization-field interaction, or by a spin-polarized 
electric current 0,0,11, SEE EH, El through the so-called 
spin-transfer torque (STT), or by a laser light [Hj]. In 
terms of applications, controlled manipulation of a mag- 
netization by a magnetic field and/or a spin-polarized 
electric current is and would be dominating techniques 
in information storage industry. The magnetization dy- 
namics of a Stoner particle under the influence of a mag- 
netic field and/or a spin-polarized current is described 
by the so-called nonlinear Landau-Lifshitz-Gilbert (LLG) 
equation that does not have analytical solutions. Most 
of the results obtained so far on magnetization rever- 
sal are either from a numerical solution Q of a g iven 
field/current pulse or based on a specific feature 0,E3, U| 
of the LLG equation. Important issues are to lower criti- 
cal fields/currents required to reverse a magnetization[llj| 
and to design a field/current pulse such that the magne- 
tization can be switched from one state to another as 
quickly as possible 0, El- 
Many reversal schemes ji, 111 III have been proposed 
and examined. To evaluate different reversal schemes, it 
is important to know the theoretical limits of the criti- 
cal switching field/current, and the optimal field/current 
pulse for the fastest magnetization reversal. The answers 
to the above questions were only known recently for the 
simplest Stoner particles of uniaxial magnetic anisotropy 
when either a magnetic field or a spin-polarized electric 
current H, El is used to reverse the magnetization. Un- 
fortunately, those ideas and approaches, based on the ro- 



tational symmetry of the system around its easy-axes, are 
not applicable to a non-uniaxial Stoner particle. In fact, 
the method does not even work on an uniaxial Stoner 
particle reversed by a spin-polarized electric current and 
a static magnetic field non-collinear to its easy-axis, a 
possible useful model in current-controlled MRAM (mag- 
netic random access memory) devices. The solutions (an- 
swers) of these questions for an arbitrary Stoner parti- 
cle under combined influences of a magnetic field and 
an electric current are the themes of current work. In 
this letter, the Euler equation of the optimal reversal 
trajectory is derived. For the special case of reversal of 
uniaxial Stoner particles by either a magnetic field or a 
Slonczewski-type of STT, analytical solution of the equa- 
tion can be obtained, which accord with the early results 
in References 5 and 12. In general, the equation can be 
numerically solved easily for a given Stoner particle of 
a known magnetic anisotropy energy. Given a reversal 
trajectory, finding the required field/current pulse is an 
algebraic problem. The generality of the current theory 
is demonstrated on a biaxial Stoner particle. 

Consider a Stoner particle of a magnetization M = 
Mm under the influence of an external magnetic field H 
as well as a polarized electric current of /, where M is 
the saturated magnetization of the particle and rh is the 
unit vector of M. Theoretical studies^, 0, Q show that 
the STT r is proportional to the current with following 
form 

r d(MV), yfU , _ , ,^ 

T = — - — \stt = g{P,m- s)m x (m x s), 

at noe 

where s is the polarization direction of the current. In the 
expression, V, h, /Uo(= 47rx 10~ 7 N/A 2 ), and e denote the 
volume of the magnetic nano-structure, the Planck con- 
stant, the vacuum magnetic permeability, and the elec- 
tron charge, respectively. 7 = 2.21x 10 5 (r ad/ s) / (A/ m) is 
the gyromagnetic ratio. The exact microscopic formula- 
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tion of STT is still a debating subject @,@|- Most of exist- 
ing theories differ themselves in different ^-functions that 
depend on the degree of the polarization P of the current 
and relative angle between m and s. Many experimental 
investigations [9( so far are consistent with the result of 
SlonczewskiQ, g = 4P 3 / 2 /[(l + P) 3 (3 + m ■ a) - 16P 3 / 2 ], 
which will be assumed in this study whenever we need 
an explicit expression of g. 

The dynamics of the magnetization M is governed by 
the generalized LLG equation 12| , 



dM - dM , , , 

— — = -jM x Ht+amx —— + jaiM x (M x s), (1) 

where a is the phenomenological dimcnsionlcss damp- 
ing constant and aj = hi g / (jiQeM 2 V) is a dimensionless 
parameter of Slonczewski STT^]. H t is the usual total 
effective magnetic field. Notice that M is a constant ac- 
cording to Eq. (jTJ), it is convenient to re- write the LLG 
equation in the dimensionless form, 
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dm 
] ~dt 



-m x hi — rh x (m x h 2 ), (2) 



where hi = ht + aais, and h 2 = ah t — ajs. t in 
Eq. |(3J) is in the units of (■yM)~ l . Both magnetiza- 
tion and magnetic field are in the units of M. The to- 
tal field h t = Ht/M = h + hi includes both the ap- 
plied magnetic field h — H/M and the internal field hi 
due to the magnetic anisotropic energy density w(fh), 
hi = —Wifiw{m)/^Q. hi and hi are in general non- 
collinear, thus the dynamics with the additional STT 
term in Eq. ([T]) is quite different [12j from that without 
this term which describes a Stoner particle in a magnetic 
field only: The particle energy can only decrease in a 
static magnetic field since the field cannot be an energy 
source 14). However, a polarized electric current can pump 
energy into a Stoner particle through the STT. Thus, 
STT allows even a dc current to be an energy source, 
and dynamics and physics of a Stoner particle under a 
arc much richer than that under a static 
magnetic field. According to Eq. the magnetization 
undergoes a precessional motion around field hi and a 
damping motion toward field hi. 

In terms of polar angle 9 and azimuthal angle </> of m 
in the spherical coordinates, Eq. ([2]) becomes 

a )9 = aj(as,p — sg) — a(- 
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(1 + a 2 )sin6»< 
dw 
~d9 
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h e ) + 



1 dw 



sin t 



Fi, 



aj(as 6 + s^) - he + 
1 dw, 



(3) 



sin0 dej) 



)=F 2 



where / means derivative of / with respect to time t. 
Here s r , sg, s</> are the r, 6, and <f> components of s, and 



Since the radial component of the external magnetic field 
does not appear in the dynamics equation of the magne- 
tization. h r shall not affect the magnetization dynam- 
ics, and we shall always assume h r — for an optimal 
field pulse. However, for STT induced reversal, s r will 
affect the magnetization dynamics through the STT co- 
efficient ai which is a function of P and s r . In general, 
Fi, (i — 1,2) are functions of 6, 4>, he, h$, se, s$, 
and s r - The switching problem is as follows: Before ap- 
plying a field/current, magnetization rh has two stable 
states, mo ('north pole' 9 = 0) and —mo ('south pole' 
9 = tt) along its easy axis (z-axis). Initially, the particle 
is in state m , and the goal is to use a field pulse and a 
spin-polarized electric current pulse to switch the mag- 
netization to — m quickly. Both the field direction and 
the spin polarization direction may vary with time. 

It will be beneficial to make a qualitative description of 
the relationships among field/current pulses, reversal tra- 
jectory, and reversal time. Given a field/current pulse, 
magnetization m will move following Eq. ^ if m is 
no longer a stable state under the field/current pulse. 
Eventually, the system will end up at one of its new sta- 
ble states[3j. If the pulse is proper, its evolution path 
4>{9) may pass through —mo. When that happens, the 
path is called a reversal trajectory, and the corresponding 
field/current pulses are reversal pulses. Given a reversal 
field/current pulse, there is a unique reversal trajectory 
which is the solution of Eq. ([2]) , a nonlinear problem with 
no general analytic solution. However, knowing a reversal 
trajectory (f)(9), which pass through 9 = 0, tt, there are 
infinite number of reversal pulses satisfying, according to 
Eq. ©, 



6' sin 9Fi - F 2 = G 2 = 0, 



(5) 



where <fi' = d(f>/d9. All hi and Sj (i — r,9,<p) satisfying 
Eq. ([5]) form reversal pulses. Thus, there is a functional 
relationship between reversal pulse and reversal trajec- 
tory. What is more, given a reversal trajectory, magneti- 
zation "velocities" 9 and (f> are linear in h and I. In other 
words, the reversal time shall be monotonic in / and h, 
magnitude of h. As h and/or / approach infinity, the 
reversal time goes to zero. Thus, the optimization prob- 
lem (shortest reversal time) is meaningful only when one 
imposes additional conditions on h and /. One sensible 
way is to fix / and h, i.e., 



Go 



0. 



(6) 



According to Eq. ([3]), the reversal time T along a re- 
versal trajectory is given by 



T 



d6 



* (1 + a 2 )d9 



(7) 



1 = 0. 



(4) 



The optimization problem here is to find the optimal re- 
versal trajectory for the fastest magnetization reversal 
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under the constrains of Eqs. |@J and ([6]). Thus T is 
minimal against the variation of the reversal trajecto- 
ries (f)(9) and field/current pulses hi and Sj (i — r, 6*, </>) 
that are linked together by Eq. ([5]). This is just the 
Euler-Lagrange problem. Using Lagrange multipliers 
method, the optimal reversal trajectory and optimal re- 
versal field/current pulses are given by 

dF 

<nr ! " 

dF 

dhi ' 



d dF 
d dF 
d dF _ dF 



(i = e. 

(i = r, 



(8) 



F 



1 



F] 



XiG\ + X 2 G 2 + A3G3, 



Gx = 0, G 2 = 0, G 3 = 0, 

where ' — jg, thus, /' means derivative of / with respect 
to 9. Ai 2,3 are the Lagrange multipliers that can be 
determined self-consistently from Eqs. (|8j). 

Eqs. ([8]) are the central result of the current work. For 
a given /, P, h, and particle (w(rh) is known), the solu- 
tion of the equations yield the optimal reversal trajectory 
4>(9) and optimal reversal field/current pulse hi and Sj 
(i = r, 9, (j)). Knowing the information, it is a straight- 
forward task to find the theoretical limit of minimum 
switching field/current. Obviously, one can simultane- 
ously optimize the field and the current in a controlled 
manipulation. From an application point of view, such a 
set up could be complicate and expensive, and one may 
be more interested in using one means only, i.e. either a 
field or a current, to manipulate a magnetization. Two 
cases may be interesting. Case A): Fix / and s, and vary 
hff and h^. Reference 5 is a special example of the case 
with 1 = 0. Case B): Fix the magnetic field h and the 
magnitude of current I, and vary the current polariza- 
tion. Reference 12 is a special example (h = 0) of the 
case. 

In order to demonstrate that Eq. |(5J) is capable to find 
the optimal reversal trajectory for an arbitrary Stoner 
particle, let us consider case A), magnetic field induced 
magnetization reversal, aj = 0. After some tedious and 
straightforward calculation, the Euler equations become 



d 2 L 



d 2 L 



dw a 
h ~~d9~sint 



d 2 L 



dw 



j. dw 

J2 = -a^rrr 



1 dw 



d9 sin 9 d<j) ' 



L = 



± 



fl + /I - (1 + o^W 

J{1 +^' 2 sin 2 0)(1 + a 2 )h 2 - (/1 - (f>> sin9f 2 ) : 



(9) 



where /i ( 2 are functions of 9 and <f> determined by mag- 
netic anisotropy energy w(m). w = —km 2 /2 corresponds 



to References 5, and it can be shown that Eqs. ([9]) is 

2 cos 9 



t(h/ h c -sm20)' 



(10) 



Its solution is 



»0 -l- — 



■In 



h-h 
^Ah/h. 



arctan[ 



2h P . 



h — h r 



-)] 



+ l)/2 + cos(6» - f) 



2ay/h + h c y/(h/h c + l)/2-cos{9 - f) 

where <po is an arbitrary constant reflecting the </>— 
symmetry of the system. This is exactly the same 
optimal trajectory equation obtained in Reference 5 
with the theoretical limit of switching field being h c = 
max{ak cos 6 sin 9/y/l + a 2 } = ak/(2y/l + a 2 ). All 
other results in Reference 5 can be re-derived from this 
optimal reversal trajectory. 

For a biaxial Stoner particle of w(9, (f) = 
-(fci/2) cos 2 9 + (k 2 /2) sin 2 9 cos 2 0, Eq.® may not 
be easy to solve analytically, but numerical solutions 
are easy to find. The solutions for a = 0.1, fci = 1, 
and h = 1(3> h c , which is about 0.05 for k 2 = 0) are 
presented in Fig. 1 for fc 2 = 0,1,5,10, respectively. 
Fig. la shows the optimal trajectories. It is interesting 
to note that the optimal trajectories for k 2 = has 
rotational symmetry around z-axis (an infinite number 
of solutions). The one shown in Fig. la is just a 
particular solution with <fi(t = 0) = 0. For k 2 7^ 0, the 
rotational symmetry is broken, but there are still four 
equivalent trajectories due to the biaxial symmetry. The 
ones whose inital 0's are in the range of [0, w/2] are used 
in the figure. Fig. lb and lc are the time evolution of 9 
and (j> along the trajectories. It is surprising to note that 
switching time is shorter and shorter as k 2 increases, 
and the potential landscape is less smooth. Fig. Id and 
le are the corresponding optimal reversal pulses. 

(a) 





- k 2 =0.0 
-k ; =1.0 

- k =5.0 



k=10.0 (©) 



FIG. 1: Numerical solutions of biaxial model with a = 0.1, 
fei = 1, h = 1, and various &2 = 0, 1, 5, 10. a) The optimal 
reversal trajectories, b) 6 vs. t. c) <f> vs. t. d) he vs. t. e) h^, 
vs. t. 

Similarly, one can obtain all the results in Reference 
12 from Eqs. {SJ. Only a spin-polarized current is used 
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to manipulate the magnetization reversal of a uniax- 
ial Stoner particle there, and the system specifities are 
h = 0, I = const., P = const., and w — w(cos9) with 
—dw/d(cos9) = f(cos9). After some algebras, the opti- 
mal trajectory equation is 



/(cos5) 



aQ[I/I c - /(cos 9) wa.6/Q] ' 



(11) 



where I r 



_ ix a eM 2 V 



— =yQ is the theoretical limit of 

\/l+a z 

the critical switching current, C 



hG(P) 



= max{f (cos 9) sin 9} 
for 9 € [0, 7r] depends on the magnetic anisotropy energy 
landscape, and G(P) — g(P, s*)y/l — s* 2 comes from the 
g-function in Slonczewski STT. The corresponding opti- 
mal polarization pulse is given by 



(1 + P? 



,*2 



16P 3 / 2 - 3(1 - py 



-as a 



(12) 



The trajectory equation has exactly the same form as 
that of magnetic-field case, reflecting the fact that both 
cases are described by the same dynamical equation ((2j. 
These are exactly what were found in Reference 12 from 
a special observation. 

Our derivation of the Euler equation is based on the 
following observations. 1) Given a field/current pulse on 
a macro-spin of a magnetic nano-structure, the trajec- 
tory of its magnetization is uniquely determined by the 
corresponding LLG equation 2) Knowing a reversal 
trajectory <f>(8) (not 9(t) and cf>(t) ), there are an infinite 
number of possible pulses, satisfying condition ([5]), that 
can reverse the magnetization along the trajectory. 3) 
The issue of the optimal magnetization reversal trajec- 
tory is an Euler-Lagrange problem with constrains. 

The Euler equation is more useful than the LLG equa- 
tion in designing an optimal field/current pulse, although 
the Euler equation is derived from the later one. With 
only a LLG equation, one can only determine one mag- 
netization trajectory for one specific field/current pulse. 
Thus, one does not even know whether the pulse can re- 
verse a magnetization! Before obtaining all solutions of 
all possible pulses, one would not be able to tell which 
one is the best. With the Euler equation, one can ob- 
tain the optimal reversal trajectory and reversal pulses 
directly for a given particle. It should be pointed out 
that the current work address the same issues as those 
in References 5 and 12, but it provides a general and 
unified theory which can deal with combined effects of a 
field and a current for an arbitrary Stoner particle. 

In conclusion, a unified theory for the optimal mag- 
netization reversal trajectories and reversal field and/or 
current pulses of an arbitrary Stoner particle is presented. 
The theory provides the Euler equation of the optimal re- 
versal trajectory along which the magnetization reversal 



is the fastest is derived. Our early results on the critical 
switching field/current for a unixial Stoner particle are 
reproduced with the unified theory. The optimal magne- 
tization reversal of a biaxial Stoner particle is also solved. 
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